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Abstract 

We investigate test particle trajectories in warped spacetimes with a thick brane warp factor, a 
cosmological on~brane hne element and a time dependent extra dimension. The geodesic equations 
are reduced to a first order autonomous dynamical system. Using analytical methods, we arrive 
at some useful general conclusions regarding possible trajectories. Oscillatory motion, suggesting 
confinement about the location of the thick brane, arises for a growing warp factor. On the other 
hand, we find runaway trajectories (exponential-like) for a decaying warp factor. Variations of the 
extra dimensional scale factor yield certain quantitative differences. Results obtained from explicit 
numerical evaluations match well with the qualitative conclusions obtained from the dynamical 
systems analysis. 
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I. INTRODUCTION 



ield, including its effects, 
, [Sj associated with the 



Much of the initial interest and subsequent research on the hypothesis that we live on 
an embedded, timelike submanifold (the brane) of a higher dimensional [D > 4) Lorentzian 
spacetime (warped or unwarped) was focused on the assumption that the scale of the extra 
dimension is independent of the coordinates. The original work of Randall and Sundrum 
(RS) on warped braneworlds, published a decade ago, did refer to the idea of the scale 
of the extra dimension being spacetime dependent, while addressing the issue of stability, in 
a two-brane scenario. In such a model, it is necessary to stabilise the inter-brane distance 
(modulus) to a fixed value, in order to avoid a collapse of one brane onto another. A 
spacetime dependent bulk scalar field, was used to achieve the desired stability of the two- 
brane system, through the Goldberger-Wise mechanism A good deal has since then 
been done with the assumption of a spacetime dependent radion 
through various types of couplings, on particle phenomenology 
two-brane RS model. 

In a single brane scenario or from a purely higher dimensional bulk perspective, the 
space-time dependence of the metric function(s) associated with the extra dimensional co- 
ordinate (s) basically imply that the scale of the extra dimension depends on the on-brane 
or four dimensional spacetime coordinates. To visualise this, it is easiest to go back to the 
standard Kaluza-Klein (KK) universe models with different scale factors associated with the 
evolution of each set of non-compact /compact dimensions (usual or extra). The difference 
between today's warped braneworlds is the warped geometry and also the non-compact 
extra dimension(s). We shall be concerned with single brane scenarios in this article. 

Since the early days of General Relativity, one possible way of understanding the nature 
of the gravitational field has been to study geodesies and geodesic deviation 6]-|8] in the 
given background geometry. While the spacetime variation of the connection reflects on 
the trajectories, the effects of curvature variations clearly control geodesic deviation. This 
is true in any dimension and in any metric theory of gravity. Usually, generic statements 
(eg. existence of orbits of different types, oscillatory/exponential behaviour in trajectories, 
focusing and defocusing etc.) and their relation to variations in the metric functions, are 
difficult to extract, though the primary goal in such studies (on geodesies) must always be 
so. 
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Over the years, several authors [9|, [id, lUl have investigated the behaviour of geodesies 



in background geometries of various higher dimensional theories of gravity. Recently, in 
[l^ . it was shown briefly, that in warped product spacetimes it is possible to have classical 



12|, 



confinement of test particles. In our work here, we confirm and add to the results in 
through detailed analysis, both numerical and analytical. 
Let us assume a bulk line element of the form, 

ds^ = e^/W [-dt^ + a'it) dX^] + b\t) da\ (1.1) 

where dX^ = dx^ + dy^ + dz^ . Here, the function h{t) represents the scale of the extra dimen- 
sion while the a{t) and e^^^"^ are the usual cosmological scale and warp factors respectively. 
In what follows, we shall largely be concerned with the effect of the geometric properties 
of the bulk spacetime on the trajectories that can exist in them. We intend to delineate 
how features change when we vary the nature of each of the three functions a(t), h{t) and 
/(cr). The geodesic equations cannot be solved analytically (modulo a few simplistic cases) 
in a spacetime as comphcated as represented by (11.11) . Thus, we make use of the dynamical 
systems approach as well as numerical methods, through which we are able to extract some 
useful information regarding the behaviour of trajectories. For numerical computations of 
the differential equations involved in this work, we have used standard numerical codes, with 
appropriate initial conditions on the variables associated with null and timelike geodesic evo- 
lution. 

The paper is organised as follows. In the next Section, we elaborate on the spacetime geom- 
etry and discuss our choices for the functions a(t), h{t) and /(cr) which we use henceforth. 
Then, in the subsequent Sections, we discuss the nature of geodesies in detail for our specific 
choices of the functions. Finally, we make a comparison of our results with those in certain 
limiting scenarios and conclude with a summary of the results obtained. 



II. THE BACKGROUND SPACETIME GEOMETRIES 

We consider the metric given by equation (11.11) rewritten in conformal time (?/), as follows, 

ds^ = e^^^''^a\r])[-d7]^ + dX.^] + b^{7])da^. (2.1) 

To arrive at concrete results in the following Sections, we need to choose the functional 
forms of the warp factor, cosmological and extra dimensional scale factors. For the warp 
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factor e^-^ we choose, 

{—log (cosh ka) — > represents a decaying warp factor, 
log (cosh A;cr) represents a growing warp factor, 

which correspond to the well-known thick brane models Q, [l^. In such models, the brane 
is dynamically generated as a scalar field domain wall (soliton) in the bulk. Note the warp 
factor in such models is a smooth function of the extra dimension, unlike the RS case where 
we have /(cr) = — A;|o"| (i.e. a function with a derivative jump-such branes are called thin 
branes). We choose to work with thick branes mainly to avoid the jumps and delta functions 
which will appear in the connection and curvature for thin branes. 

What do we choose for the scale factors? Two different combinations of a{ri) and b{ri), 
chosen as models to represent different kinds of time evolution, are given below in terms of 
the cosmological time "t" , 

{a{t),b{t)} = {i) {aot•'\bo + b^t-''}, and (ti) {a[,e^(*-*«), 6[, + 6;e-^^(*-*°)}, (2.3) 

where we take ui to span over an open interval (0, 1), so that in (z) a{t) represents an 
expanding but decelerating on-brane line element (which is radiative for ui = ■^), and in (ii) 
we take H to be positive, representing an accelerating de-Sitter on-brane line element, bo, 
bi, 1^2, b[ and (3 are positive so that we have a decaying extra dimension which stabilizes 
to a finite value as t — > oo. The type of metric functions we consider here may be obtained as 
analytic solutions of the five dimensional Einstein equations with the corresponding Einstein 
tensors providing matter energy-momentum profiles through the Einstein field equations. As 
shown in [l^, for similar metric functions, the corresponding matter stress-energy can indeed 
satisfy the Weak Energy Condition. 

In the first case, the conformal time = J We choose to work with an initial 

condition such that at t = to, rj = rjo = 1 which, in turn, implies ao = and ^7 = ( ) 
where the domain of ?7 is 1 < < oo. We also assume that evolution of our universe begins 
at t = to (oi' ^ = 1) when the scales of ordinary and extra space dimensions were same (i.e. 
(^(to/rjo) = bito/rjo)). In terms of conformal time, we are led to, 

to bi / — \ 

ci'i'n) = V^'"^ and b{i]) = TU^ ('^ - V^'"^ ■ (2-4) 

i — Ui i — i^l tg 

The stabilization condition further implies, 

bo = j^-%>0 or b,<^. (2.5) 
1 - z/i to 1 - 



In order to reduce the number of free parameters, we consider bi = tg^'^^ (which is, in fact, 
the maximum possible value for bi) as z/i can only be close to zero (but not equal to zero). 
Therefore we have. 



Using a similar prescription for the de-Sitter model, we obtain. 



(2.6) 



(2.7) 



where b'^H < 1 and < ?7 < 1. Let us consider, ui = V2 = \, to = H = 1, b'l = \ and 
(3 = 1 which leads to the following two different combinations of the scale factors, 

(A) a{'n)=2r,, 6(ry) = 1 + i, 

(B) a{^) = ^, %) = 1 - f . 

These combinations of scale factors will henceforth be abbreviated as set (A) and set (B). 
Fig. [1] gives a pictorial representation of the warp factors (decaying and growing) and the 
scale factors corresponding to set (A) and set (B). 
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FIG. 1: Variations of warp factors and the scale factors for set(A) and set(B) respectively. 



III. GEODESICS 

The general form of the constraint (i.e. qabu^u^ = — e) for null and timelike geodesies 
corresponding to the metric (12.11) is given below, 

e^^^'^^a^iT]) [-7)2 + X^] + b{r]y a' + e = 0, (3.1) 

where e = 1 and denote the cases corresponding to the timelike and null geodesies respec- 
tively and a dot here represents the differentiation with respect to the afiine parameter A. 
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As Xi (such that xi = x, X2 = y, and X3 = z) is cychc, an obvious conclusion for the metric 
(iH) is, 

Xi = 3.2 

where Dj's are integration constants. The effective geodesic potentials V„ and \4 are then 
calculated as, 

,,,,,_ ^ + a\v)e'f^''K-v' + ^i') .... 



.2 e + &2(T2 



r)2 -4/(<7) 

(3.4) 



2a4(r/) 



e2/Wa2(?^)_ 

where we have considered x = y = z, which will be used throughout hereafter. In order to 
find K-(cr) (or V^(x)), we shall make use of solutions for x{X), ri{X) and a{X) in the right 
hand side of equation (13. 3p (or equation (13.41) ) and, subsequently, make a parametric plot 
of K-(A) versus cr(A) (or \4(A) versus x{X)) to arrive at a pictorial representation of Vo-(cr) 
(or V"^(x)). It may be noted that, we have absorbed the factor "e" in the effective geodesic 
potential. Therefore, for both timelike and null geodesies, we are essentially looking at test 
particles with zero "total energy" (from a mechanics point of view). 



A. Dynamical systems analysis 

The geodesic equations corresponding to the metric ( 12. ip derived for general warp and 
scale factors turn out to be, 



[^^ + E,,(i.)']+e-^^('^) ^^a2 + 2/'(a)7)a = 0, (3.5) 



2a.'(??) , , V 

Xi + —^^fjx, + 2f'{a)x,& = 0, 3.6 
a{ri) 

here the dots represent the differentiation with respect to the affine parameter A (or an 
arbitrary parameter, for null geodesies) while primes denote differentiation of the functions 
with respect to their corresponding independent variables, r] or cr. The geodesic equations 
corresponding to x, y, and z can be reproduced from equation (13. 6p and they all have an 
identical structure except for the constants Di. The full geodesic equations in a spacetime 
corresponding to the metric (12. ip are very difficult to solve analytically. However, the Eqs 
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(13. 5p - (13. 71) can be recast as the following dynamical system of first order, coupled differential 
equations, 

f] 

Xi 
& 

and X 

It will be interesting to see if one can analytically draw some useful conclusions about the 
behaviour of the geodesies by using the above set of equations. For the growing warp factor, 
/'((j) ~ tanh(/i;cr) and if we expand this function about the location of the brane (i.e. cr = 0), 
the dominant term in the neighbourhood of that point will be linear in a. In such a scenario, 
for the timelike case (e = 1), the Eqs (I3.10p and (13. lip resemble the following equations, 

c"" ~ X ^ positive terms, and x~ — o" x positive terms, (3-12) 

which are qualitatively similar to a dynamical system representing simple harmonic motion 
(SHM) {X = P, P = —X} (modulo a multiplicative factor). From this consideration, one 
can intuitively argue that the Eqs (]3.8p - (]3.1ip will have oscillatory solutions for a and x 
(though, not simple harmonic). On the other hand, for a decaying warp factor, /(cr) = 
— log(cosh ka), the above set of equations given by (I3.12p lead to a system with exponential 
solutions. Thus, there should not be any oscillatory trajectories with decaying warp factors. 

From equation (I3.12p . it is also evident that, for any growing warp factor one does not 
have oscillatory timelike trajectories. For example, in the Randall-Sundrum scenario, we 
have, f{a) ~ |cr| and /'(cr) is therefore a step function, which cannot have an expansion 
with a term linear in cr. So there exist no oscillatory solution for this case. 

In order to further analyse the behavior of trajectories analytically, let us look at the 
specific case where 6(77) is a constant (say 6(77) = 1). Since we are particularly interested 
about motion in the bulk, we choose, Xj = (as Xj does not appear in the other three 
equations i.e. in (13. Sp . (I3.10p . and (13. lip . With these considerations, we arrive at following 



a{r]) y IP'ir]) 
X 



-/'(^) ( e 



X^ 



(3.8) 

(3.9) 
(3.10) 

(3.11) 
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set of equations, 



V = -^V^ + X^ (3.13) 
a{r]) 

a = X, (3.14) 

and X = -f'{cT){e + x'). (3.15) 

For timelike geodesies and growing warp factors, from Eqs (13.141) and (13.151) . we obtain, 

~r = ~7^, — — T~ =^ Cisecha = ^/lTx^, (3.16) 
dx (1 + X ) tanner 

where Ci is an integration constant. Similarly, Eqs (I3.13p . (13.151) and (I3.16P further lead to, 
^ = -2v ^Jcf-l-x' =^ x=\lcl-U\n{C,~r,^), (3.17) 

where C2 is another integration constant. Now combining the Eqs (13.161) and (I3.17p . one 
easily obtains, 

a = sech~^ 



1 + (C?-1) sm'{C,-^^)\. (3.18) 

The argument of the "sech~^" function contains a sinusoidal part which guarantees an 
oscillatory behaviour of cr(?7), though its details depend on the parameters (Ci and C2) 
involved in the analysis. It is worth mentioning here that even if we let airf) to be a 
constant, oscillations persist (only the factor sin^(C2 — rf) gets replaced by sin^(C2 — 2?])). 
This confirms that the oscillations are solely due to the presence of a growing warp factor. It 
has been checked that, for a decaying warp factor, no such analytic expression can be found. 
In fact, the Eqs (13. 8p . (I3.10p and (13.111) form a dynamical system of their own, irrespective 
of equation (13. 9p . Unfortunately, with e = 1, this system does not have any fixed points for 
a finite value of ri. But, for a static extra dimension, we can investigate the solution space 

n 

[16l | for the following set of two equations, 

a = Xi and X = ±(1 + x^) tanh fccr, (3.19) 

where + and — signatures denote the decaying and growing warp factors respectively. It 
can easily be shown that the point (a, x) = (0, 0) is a saddle point in the phase space for 
the decaying warp factor and a centre in case of a growing warp factor. Thus, we can 
expect oscillatory solutions for cr(A) and x(A) in the latter case (i.e. for a growing warp 
factor). We present the phase space plots of a and x with two different warp factors in 
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FIG. 2: The solid lines indicate different flows in the solution space of for differ- 

ent initial conditions, such as, in the left plot (for decaying warp factor): {cj(0),x(0)} = 
(0.25,0), (0,-0.25), (-0.25,0), (0,0.25) (here A runs from -1.3 to 1.3), and in the right plot (for 
growing warp factor): {cj(0), x(0)} = (0, 1), (1, 2), (2, 3) (here A runs from to 6.5). 

the Fig. [21 A decaying warp factor only helps a particle to move away from the brane, 
except for the initial conditions satisfying a -|- x = (where the particle tends toward the 
location of the brane with an ever-decreasing velocity). This essentially means that only 
for the abovementioned initial conditions, particles will end up accumulating near the brane 
with increasing affine parameter. Particles with initial condition a — x = will be repelled 
away from the brane. There are two other kinds of trajectories. The trajectories in the 
quadrants (a -f x > 0, cr — % < 0) and (a -f x < 0? — x > 0) may cross the brane at some 
point of time (depending on the initial conditions), and the trajectories in the quadrants 
(cr-|-x>0,cr — x>0) and (a-|-x<0,cr — x<0) will bypass the location of the brane 
(though they may come very close, depending on the initial conditions). On the other hand, 
the closed curves, in the case with a growing warp factor, indicate that, cr(A) and xi}^ have 
an oscillatory behaviour. One can also notice that the nature of the oscillations depend 
heavily on the initial conditions. A particle, initially on or close to the brane and having a 
small velocity, executes an almost sinusoidal oscillation. This oscillatory pattern will tend 
towards a square wave profile as the distance of the launching point from the brane, or the 
initial velocity, increases. The phase portraits in Fig. [2] are for a static extra dimension. 
We expect to see a similar solution space for a dynamic extra dimension {b{r])) too when 
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1] ^ oo because according to our choices 6(77) gets stablised to a constant value for large rj. 
The closed trajectories for a static extra dimension and a growing warp factor would become 
limit cycles in phase space when we have a dynamic extra dimension. 

One may also note that for null geodesies, the eigenvalues of the linearised system at the 
origin in the phase space vanishes. It is therefore not possible to comment on the solution 
space for this case. 

B. Numerical evaluations 

The full geodesic Eqs ( I3.5p - fl3.7l) for the null and timelike geodesies can obviously be 
solved numerically. We present some features of these geodesies, graphically, with the set 
(A) and (B) for the scale factors and with the growing and decaying warp factors of the 
thick brane model. All of these cases are presented in Figs. [3] -|H] where the figure captions 
are self-explanatory and describe the corresponding results. We have shown the solutions 
for timelike and null geodesies in the case with a growing warp factor. Here, we have 
observed remarkable differences between the null and timelike cases (see Figs. IH [5l [7] and 
[8]). Otherwise, where the differences between the results for timelike and null cases are small 
(i.e. results remain qualitatively similar for the case of a decaying warp factor), we present 
only the timelike geodesies (see Figs. [3] and El). The initial conditions for our numerical 
evaluations are chosen asi; = ?/ = i = (T = 0.1atA = and ?) at A = is calculated by 
using the null and timelike constraints which differs between models. 

We have observed that for a growing warp factor, the timelike geodesies have an oscillatory 
behaviour for both the cases (i.e. set (A) and (B)) as presented in Figs. [Hand [7] respectively. 
Such an oscillatory behaviour is absent in the trajectories of null geodesies in all the cases. 
This essentially means that massless particles are free to escape into the bulk. The oscillatory 
behaviour of the timelike geodesies for a growing warp factor also manifests itself in terms of 
the potential V{a) as shown in Figs. [Hand [7] for set (A) and set (B) respectively. The fact 
which is worth mentioning about these trajectories is that a massive particle never moves 
away to an infinite distance from the location of the brane. It seems that the trajectories are 
automatically localised on or near the brane in the presence of a growing warp factor. The 
potential in this case is in fact that of an anharmonic oscillator where both the amplitude and 
frequency of the oscillations are monotonically increasing and converging towards a limiting 
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FIG. 3: Timelike geodesies for set (A) with a decaying warp factor. 



5 10 15 20 

A 



v[x) 


-0.001 




0.05 0.1^ 


J^r^ DT2 


0.25 


-0.002 










-0.003 










-0 . 004 










-0.005 












FIG. 4: Timelike geodesies for set (A) with a growing warp factor. 



value with an increasing A. We will again return to the interpretation of this feature in terms 
of the warp factor and the scale of the extra dimension, in the Section IIVI Our numerical 
findings in this section tally well with the results found earlier in Section III, via the analysis 
of an autonomous dynamical system derived from geodesic equations. 
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FIG. 5: Null geodesies for set (A) with a growing warp factor. 
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FIG. 6: Timelike geodesies for set (B) with a decaying warp factor. 

IV. EFFECTS OF WARPING AND TIME DEPENDENCE IN EXTRA DIMEN- 
SION SCALE AND CONFINEMENT OF TEST PARTICLES 



In order to figure out the effects of warping and a time dependent extra dimension, we 
will compare the generalised braneworld scenario (discussed in the previous section) with 
the situations where we do not have any warp factor or a time dependence in the extra 
dimension scale. Let us first consider, the metric (12.11) without the warping factor. This 
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FIG. 7: Timelike geodesies for set (B) with a growing warp factor. 




FIG. 8: Null geodesies for set (B) with a growing warp factor. 

leads to the following constraint, 

a^T]) [-drf + rfX^] + h{r]fa^ + e = 0. (4.1) 

Now cr becomes cyclic in the geodesic Lagrangian and gives us Vx{x) = —D'^/{2a^) and 
V„{a) = — P^/(26'^) which eventually lead to, 




With a nonzero P, we could not find an analytic solution for Eq. 14.21 (though it is exactly 
solvable for P = 0). We have therefore solved the geodesic equations for P 7^ numerically. 
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In the case where b{T]) is constant, the constraint (13.11) gives, 

The Eq. (14.31) is also not analytically solvable except for the case with D = 0. As above, we 
solve the geodesic equations numerically for the general situation when D ^ 0. 

In order to see the specific differences between the general case (eg. timelike geodesies for 
set (A) with growing warp factor as shown in Fig. H]), and two of its subcases (one without 
a warping factor and the other with a static extra dimension), the plots of the potential 
Va{cr) and the trajectories for both of these two subcases are presented in Fig. [9l Here, 
we have considered D = 1 and P = 1 for the purpose of numerical computations. The 
differences between the two subcases are quite clearly visible. We have also worked through 
other similar cases, which we do not mention here. It is clear that the oscillatory behaviour 




FIG. 9: Plots in the upper and lower rows show the timelike geodesies for set (A) (with b{r]) = 
constant) for models without any warping factor and with growing warp factor respectively. 

is due to the presence of a growing warp factor (as already argued in Section III). If we 
compare Fig. M with Fig. HI where scale of the extra dimension is dynamic, the increase 
in the amplitude and frequency of the oscillations of a, may be understood as an effect of 
the decay in the scale of the extra dimension b{ri). In fact, the scale 6(77), as mentioned 
before, can be viewed as the length of an anharmonic oscillator which tends to become a 
harmonic one with an increasing value of the affine parameter. In other words, this problem 
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is equivalent to a problem of a particle in a one dimensional box whose size is decreasing 
and tending to a finite value. This is exactly how one can understand the increase in the 
amplitude of the oscillation in o"(A), in Fig. IH and its eventual stabilisation to a limiting 
value. Such features are also evident from the behaviour of the corresponding potentials as 
given in the Figs. H] and M 

Let us try to understand the above-mentioned confinement from a different point of view. 
It was shown in [12;], that in the case of a five dimensional bulk with a static extra dimension 
and a growing warp factor (which is indeed a solution of five dimensional Einstein equations 
in presence of a negative bulk cosmological constant), we have confinement of test particles. 
Further, in |l2|, it was shown that confinement will be achieved if the following energy 
condition holds, 

(°i?ABef ef - °i?,/3)M"M^ > 0, = 0,1,2,3,4; a, /5 = 0,1, 2, 3 (4.4) 

where, = are the four basis vectors and {x^} and {y"} are five and four dimensional 
coordinate systems respectively. Rab and RajS are Ricci tensors in the bulk and on the a = 
const, hypersurface respectively. The above condition implies that the local gravitational 
density of five dimensional bulk matter, as measured by an observer freely falling along a 
hypersurface o" = (or the brane), has to be greater than the gravitational density of the 
effective four dimensional matter. For the five dimensional geometry given by Eq. 12.11 the 
above constraint leads to 

(4.5) 




It may be noted that for a static extra dimension, i.e. when 6(77) ~ constant, the following 
condition, 

/" + 4/'2 > (4.6) 

is sufficient for confinement of test particles. This condition is clearly satisfied for the 
functional form of the growing warp factor we have considered. 

Now the constraint Eq. 14.51 is always satisfied for a monotonically growing warp factor 
with a growing a{ri) and a growing but decelerating b{ri). In that case, if b{ri) does not 
stabilise, then as r] 00, the depth of the potential well tends to zero (as b{ri) appears 
in the denominator of Eq. 13. 3p . As a result, confined trajectories will become more and 
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more unstable against perturbations and particles should disappear from the brane. This 
observation should rule out the models with growing and divergent b{r]), for example, some 



of the power law solutions (for b(t)) found in 15|, in presence of different bulk matter fields as 



such. On the other hand, if b{r]) grows but stabilises to a finite value a.s 1] ^ oo, trajectories 
are confined, though become less strongly bound as the depth of the potential decreases 
(towards a limiting value). 

The last two terms in the left hand side of 14. 51 are always negative for a growing a{ri) and 
a decaying b{ri). Therefore, it is obvious that the inequality 14.51 will not always hold for any 
set of metric functions (though it happens to be the case for the metric functions we choose). 
However, one can argue that as b{ri) stabilises with t] oo, constraint 14.51 converges toward 
constraint 14.61 It is thus implied that confinement will eventually be achieved during the 
course of cosmological evolution. The trajectories, in this case, will also be more strongly 
bound as the depth of the potential well increases. It is worth mentioning here that, if b{r)) 
decays to a zero value, the depth the potential 13.31 tends to infinity resulting in an ever- 
growing amplitude of oscillation which eventually diverges. This aspect may help us rule 
out models with the 6(77) decaying to zero (instead of stabilising to a finite limiting value). 

One may also find a constraint on the warp factor, for confinement, in a static bulk in 
the following simple way. For the test particles (moving along the extra dimension) to be 
confined somewhere in the bulk, their geodesic potential should have a minimum at that 
point. Now the corresponding geodesic potential is given by 

Via) = (4.7) 

where C is, a positive constant, dependent on the initial conditions. V{cr) should have a 
minima at a = ctq, where /'(co) = and e~'^^^°'"^ 7^ 0, if the following condition holds 

/"((To) > 0, (4.8) 

which means the f{a) must be a monotonically growing function of a about the location 
where the potential has a minimum. For f{(r) = log(cosha), V{(r) = The 
minimum value of this potential for timelike geodesies is — (where C > 1) at cr = and 
it reaches its maximum value +| as a ^ ±00. Therefore the zero energy trajectories are 
bounded within the a values which can be obtained by equating V{a) = 0. This leads to 

/(a)=logv^. (4.9) 
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We can therefore explain the confinement of massive particles in the context of Figs H] and 
[6l On the other hand, for null geodesies, the zero energy trajectories always have runaway 
features (Figs. [S]and[H]). 

We note that whenever Eq. 14.81 is true, the Eq. 14.61 is automatically satisfied. This also 
justifies why geodesies in a thick braneworld model with growing warp factor i.e. /(a) ~ 
log(coshcr) are confined near a = 0, whereas with a decaying warp factor we don't see such 
confinement . 

V. SUMMARY AND CONCLUSIONS 

To conclude, we provide below a summary of the results obtained, in a systematic way. 

• We have studied geodesies for the warped background geometry of a thick braneworld 
model with a cosmological metric on the brane and a time dependent extra dimension, 
in detail. Our choices for the scale factors (marked as set (A) and set (B) throughout 
this article) illustrate situations corresponding to decelerating (i.e. set (A)) and accel- 
erating (i.e. set (B)) universes. The extra dimensional scale is chosen appropriately 
which stabilises to a constant value at large time. 

• The geodesic equations in such a background can be rewritten as a first order au- 
tonomous dynamical system. Analytical insights as well as some specific solutions can 
therefore be obtained. In particular, the analysis of these equations as a dynamical 
system shows the role of the warp factor (growing/ decaying) in controlling the nature 
of the solutions (oscillatory/exponential). 

• Apart from a few specific cases, the geodesic equations, in general, can only be solved 
numerically since they are highly coupled differential equations. However, we are able 
to demonstrate some qualitative features of the dynamical system through numerical 
investigations. Detailed features of the null and timelike geodesies are brought out 
for different cases with growing and decaying warp factors. It may be noted that a 
numerical estimation is always possible for any other combination of scale factors/warp 
factor. 

• We have tried to separate out the individual effects of the warp factor and the extra 
dimensional scale factor. It is found that the emergence of an oscillatory behaviour 

17 



in timelike geodesies (absent in the case of null geodesies) is essentially due to the 
presence of a growing warp factor. It seems that massive particles are constrained 
to be confined near the brane in this particular scenario (growing warp factor), while 
massless particles like photons or gravitons can indeed access the bulk. This observa- 
tion is important in the context of localisation of fields on the brane. The behaviour 
of localisation of test particles in a geometry with a growing warp factor is reminis- 
cent of fermion localisation with a growing warp factor. Even in the presence of a 
decaying warp factor, particles can also be localised though very special initial con- 
ditions are required. The effect of a dynamical (time dependent) extra dimension on 
the trajectories is found to be largely quantitative. 

• We have clearly demonstrated that the analysis of geodesic motion obtained through 
a dynamical systems analysis matches quite well with the numerical computations. 



• We have briefly examined a confinement condition obtained in [l^ for our cases here. 
We also provide our comments on aspects of this condition and show that the issue 
of confinement puts important constraints on possible dynamic behaviour of the extra 
dimension. 

We mention that it would be useful to investigate geodesic deviation and the kinematics of 
geodesic flows in a bulk geometry with a thick brane. In the latter case, we need to solve the 
geodesic and Raychaudhuri equations simultaneously, in order to obtain the evolution and 
initial condition dependencies of the kinematical quantities: expansion, rotation and shear 
(or ESR). We believe that such work may improve our understanding about the nature 
of warped extra dimensions. We hope to report on related work on the above-mentioned 
aspects, at a later stage. 
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